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We present a combined high-resolution x-ray diffraction and x-ray resonant magnetic scattering (XRMS) 
study of as-grown BaFe2As2. The structural/magnetic transitions must be described as a two-step process. At 
Ts = 134.5 K we observe the onset of a second-order structural transition from the high-temperature paramag- 
netic tetragonal structure to a paramagnetic orthorhombic phase, followed by a discontinuous step in the struc- 
tural order parameter that is coincident with a first-order antiferromagnetic (AFM) transition at Tn = 133.75 
K. These data, together with detailed high-resolution x-ray studies of the structural transition in lightly doped 
Ba(Fei_ x Co a; )2As2 and Ba(Fei-^Rha;)2As2 compounds, show that the structural and AFM transitions do, in 
fact, occur at slightly different temperatures in the parent BaFe2 As2 compound, and evolve towards split second- 
order transitions as the doping concentration is increased. We estimate the composition for the tricritical point 
for Co-doping and employ a mean-field approach to show that our measurements can be explained by the inclu- 
sion of an anharmonic term in the elastic free energy and magneto-elastic coupling in the form of an emergent 
Ising-nematic degree of freedom. 

PACS numbers: 74.70.Xa, 75.25.+Z, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

Although most of the excitement surrounding the 
discovery^ 2 - of superconductivity (SC) in the iron pnictide 
compounds has focused on the underlying mechanism, these 
compounds also provide us with a fascinating opportunity to 
study subtle, and not so subtle, aspects of the interactions 
between SC, structure and magnetism. A great deal of this 
research has focused on the 7?FeAsO (R = rare earth), or 
7777 compounds, and the As2 (AE = alkaline earth), or 

722 family. The parent 7777 and 722 systems undergo tran- 
sitions from a high-temperature paramagnetic (PM) tetrago- 
nal phase to a low-temperature antiferromagnetic (AFM) or- 
thorhombic structure. 3-9 In the 7777 parent compounds, the 
structural transition occurs at somewhat higher temperature 
than the AFM ordering. 3 In the 722 parent compounds, these 
transitions appear coincident in temperature and coupled^r— 

The consequences of electron/hole doping for the 
BaFe2As2 compound have been studied widely. For exam- 
ple, upon electron doping through substitution on the Fe site 
by Co, Ni, Rh, Pt or Pd, the structural and AFM transi- 
tions are suppressed and split, with the structural transition 
found at somewhat higher temperature. 10-18 Superconductiv- 
ity emerges over a small but finite range in dopant concentra- 
tion. In contrast, the structural and AFM transitions remain 
coincident in superconducting samples produced by K doping 
on the Ba site, P doping on the As site, and Ru doping on the 
Fe siteM^— In Co-, Rh-, Ni-, and Ru-doped BaFe2As2, neu- 
tron and x-ray measurements revealed a suppression in both 
the magnetic order parameter and the orthorhombic distortion 
below the superconducting transition temperature, which in- 
dicates the coexistence and competition between magnetism 
and superconductivity in these systems ^ 18 ' 22- — 

The nature of the structural and AFM transitions has, it- 



self, been the subject of intense scrutiny. For CaFe2As2 and 
SrFe2As2, several neutron, x-ray, /iSR and NMR measure- 
ments have reported that the structural/magnetic transitions 
are discontinuous (first-order) and hysteretic^ 8 ^ However, 
in BaFe2As2, there has been significant debate concerning 
the character of the structural and magnetic order parameters. 
Early neutron measurements on polycrystalline samples found 
a second-order magnetic transition. 9 Neutron diffraction mea- 
surements by some groups, however, described a first-order 
magnetic transition but failed to observe any hysteresis, 26 
while neutron and NMR measurements by others found a first- 
order magnetic transition and observed a large hysteresis upon 
cooling and warming. 27 ' 28 Wilson et al. initially reported that 
the magnetic and structural transitions in BaFe2 As2 were con- 
tinuous in nature, and could be described by a simple power- 
law dependence with an critical exponent consistent with a 
two-dimensional Ising model. 29 ' 30 Later heat capacity and x- 
ray work by this group on annealed samples of BaFe2As2 
found that the orthorhombic distortion appeared first as a 
second-order transition interrupted, at slightly lower temper- 
ature, by a first-order transition to the low temperature or- 
thorhombic phase. 31 As noted by these authors, this complex 
structural transition and its relationship to the concomitant 
AFM ordering calls for further investigations. 

To shed light on the nature of structural and magnetic transi- 
tions, we present a combined high-resolution x-ray diffraction 
and x-ray resonant magnetic scattering (XRMS) study of as- 
grown BaFe2As2. In Section II we provide details of the scat- 
tering experiments on the parent BaFe2 As2 compound and the 
Co- and Rh-doped compounds. In Sections III, we report the 
results of high-resolution x-ray diffraction and XRMS mea- 
surements on the parent BaFe2As2 compound and show that 
there is a small, but distinct, difference in the temperatures of 
the structural and AFM transitions. We also studied the doped 
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compounds to substantiate our conclusion that the structural 
transition is continuous in nature, whereas the AFM transi- 
tion changes from a first-order transition at low doping to a 
second-order transition for higher doping levels as the system 
passes through a tricritical point. In Section IV, we estimate 
the position of this tricritical point for the Ba(Fei_a;Co x )2 As2 
phase diagram and employ a mean-field approach to show that 
our measurements can be explained by the inclusion of an an- 
harmonic term in the elastic free energy and magneto-elastic 
coupling in the form of an emergent Ising-nematic degree of 
freedom. In Section V, we summarize our results. 



II. EXPERIMENTAL DETAILS 

Single crystals of BaFe2As2, Ba(Fei_ x Co x )2As2 and 
Ba(Fei_ x Rh2;)2As2 compounds were produced using the self- 
flux solution growth method described elsewhere. 11 Energy 
dispersive spectroscopy (EDS) was performed to confirm the 
absence of foreign elements, and wavelength dispersive spec- 
troscopy (WDS) was employed to determine the compositions 
of the Co- and Rh-doped compounds at several points on each 
sample, providing a combined statistical and systematic er- 
ror of less than 5% of the relative elemental concentration 
[e.g. 0.018±0.001 for the Ba(Fe .982Coo.oi8)2As 2 sample]. 
Temperature-dependent ac electrical resistance data (f = 16 
Hz, 7=3 mA) were collected using a Quantum Design Mag- 
netic Properties Measurement System (MPMS) with a Lin- 
ear Research LR700 resistance bridge. Electrical contact was 
made to the sample using Epotek H20E silver epoxy to attach 
Pt wires in a four-probe configuration. 

Temperature-dependent, high-resolution, single-crystal x- 
ray diffraction measurements were performed on a four-circle 
diffractometer using Cu K a \ radiation from a rotating anode 
x-ray source, selected by a germanium (1 1 1) monochro- 
mator. For these measurements, the plate-like single crys- 
tals with typical dimensions of 3 x 3 x 0.5 mm 3 were at- 
tached to a flat copper sample holder on the cold finger 
of a closed-cycle displex refrigerator. The mosaicities of 
the BaFe2As2, Ba(Fei_ 2 ,Co :E )2As2 and Ba(Fei_ 2 ,Rh x )2As2 
single crystals were all less than 0.02° full-width-at-half- 
maximum (FWHM) as measured by the rocking curves of the 
(11 10) reflection at room temperature. The diffraction data 
were obtained as a function of temperature between room tem- 
perature and 8 K, the base temperature of the refrigerator. 

To correlate the evolution of the structure with the oc- 
currence of magnetic order, both conventional x-ray diffrac- 
tion and XRMS measurements were performed on the 6ID- 
B beamline at the Advanced Photon Source (APS) using the 
same as-grown BaFe2As2 single crystal studied with the lab- 
oratory source. The single crystal was attached to a flat cop- 
per sample holder on the cold finger of a closed-cycle dis- 
plex refrigerator with the tetragonal (HHL) plane coincident 
with the vertical scattering plane. The temperature, measured 
at a sensor mounted to the copper block holding the sample, 
was stable within ±0.002 K. Care was taken to ensure that 
heating effects associated with the incident x-ray beam were 
minimized by measuring charge and magnetic reflections in 



close proximity and using the appropriate incident beam at- 
tenuation. Measurements of both the charge scattering and 
the XRMS were done at the Fe If -edge (E = 7.112 keV)M 
The incident radiation was linearly polarized perpendicular 
to the vertical scattering plane (er-polarized) with a spatial 
cross section of 1.0 mm (horizontal) x 0.2 mm (vertical). In 
this configuration, dipole resonant magnetic scattering rotates 
the plane of linear polarization into the scattering plane (ir- 
polarization), while the charge scattering leaves the polariza- 
tion unchanged. Cu (2 2 0) was used as a polarization analyzer 
to suppress the charge and fluorescence background relative 
to the magnetic scattering signal by approximately a factor of 
200. 



III. RESULTS 

A. High-resolution x-ray diffraction and XRMS 
measurements of BaFe2 As2 

In Fig.QJa), we display [£ £ 0] scans through the (1 1 10)t 
charge peak, obtained using the laboratory source, for the 
parent BaFe2 As2 compound measured with temperature steps 
of 0.25 K. Above the structural transition temperature, T$ = 
134.5 K, we observe a single sharp peak consistent with the 
tetragonal structure. Upon cooling below T5, the (1 1 10)t 
charge peak continuously broadens and, then, clearly splits 
at T = 133.75 K concomitant with the abrupt appearance of 
two additional peaks at this temperature [vertical arrows in 
Fig. QI a )] bracketing the two inner peaks. Upon further cool- 
ing, the splitting of the two inner peaks evolves continuously 
as their intensities decrease, whereas the positions of the outer 
peaks change only slowly as their intensities increase. Below 
T = 133.0 K, the two inner peaks disappeared leaving only 
the outer peaks in evidence. We note that these observations 
are qualitatively consistent with similar diffraction measure- 
ments on an annealed sample of BaFe2As2 recently reported 
by Rotundu et al. 31 although the transition temperatures for 
their annealed sample were approximately 5 K higher. 

Having described the temperature evolution of the diffrac- 
tion peaks qualitatively, it is useful at this point to intro- 
duce some labeling of the corresponding phases. The high- 
temperature paramagnetic phase is denoted as Tet-PM. An- 
ticipating the results of our XRMS study, we label the or- 
thorhombic phase that evolves continuously over a very nar- 
row temperature range below T$ [corresponding to the inner 
pair of peaks in Fig. [Tf a )] as Ort-PM. We further label the or- 
thorhombic phase that abruptly appears at T = 133.75 K [cor- 
responding to the two outer bracketing peaks in Fig. [TJa)] as 
Ort-AFM. Structurally, we assume that Ort-PM and Ort-AFM 
differ only with respect to the values of their lattice constants 
and orthorhombic distortion at a given temperature. 

Figure [TJb) describes the temperature evolution of these 
phases. Upon cooling, a second-order transition from Tet-PM 
to Ort-PM occurs at T$ = 134.5 K followed by a first-order 
transition to Ort-AFM at T N = 133.75 K. There is a region of 
coexistence between Ort-AFM and Ort-PM from 133.75 K to 
133.0 K, and only the Ort-AFM phase is observed below this 
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FIG. 1 : (Color online) (a) X-ray diffraction scans, measured using 
the laboratory source, along the [£ £ 0] direction through the position 
of the tetragonal (1 1 10)t reflection for selected temperatures in the 
parent BaFe2As2 upon cooling. The lines present the fitted curves 
using a Lorentzian-squared line shape. The two-component fit to 
broadened peaks is illustrated for T = 134.25 K. The arrows denote 
the positions of peaks associated with Ort-AFM as discussed in the 
text. At this temperature, the integrated intensity of the Ort-AFM 
peaks are approximately 5% of the Ort-PM diffraction peaks, (b) The 
orthorhombic distortion as a function of temperature upon cooling 
and warming determined from fits to the (11 10) Bragg diffraction 
peak. 



temperature. Upon warming, Ort-PM appears at 133.0 K and 
coexists with Ort-AFM up to T' N = 134.0 K, where Ort-AFM 
disappears. The orthorhombic distortion associated with Ort- 
PM decreases continuously up to Ts = 134.5 K, where Tet- 
PM is recovered. We find no hysteresis in the transformations 
from Tet-PM to Ort-PM and <0.25 K hysteresis associated 
with the appearance/disappearance of Ort-AFM. 

In order to investigate the relationship between the struc- 
tural transition and AFM ordering in this system we have per- 
formed a combined study using high-resolution x-ray diffrac- 
tion and XRMS measurements. These simultaneous measure- 
ments eliminate concerns regarding disparities in the temper- 
ature calibration of sensors for different experiments. Using 
the configuration at the APS described in the last section we 
measured the scattering at both the charge and magnetic Bragg 
peak positions for several temperatures close to the structural 
transition. In Fig. [2 a) we show a [£ £ 0] scan through the (11 
8)t charge Bragg peak at T = 137 K, well above the structural 
transition temperature. At T = 130 K, below Tm and Ts, two 
well-separated peaks were observed [Fig. ISJb)] . These are the 
(2 8)o and (0 2 8)o charge Bragg peaks of the orthorhom- 
bic phase. The difference in intensity arises from different 
populations of the domains within the illuminated volume of 
the sample. At this same temperature, a single peak is found 
at the (10 7)o, magnetic peak position for the orthorhombic 
phase, [Fig. |2fc)] in agreement with previous measurements 
of a magnetic propagations vector given by Q m = (10 l)o [(^ 
i 1)t] with lattice constants ao > &o i 9 ' 26 i 27 i 29 i 32 For simplic- 
ity, we will henceforth label all peaks with tetragonal indices. 
Therefore, the (10 7)o magnetic Bragg peak will be referred 
to as (| i 7)t, keeping in mind that the magnetic peaks are 
displaced from £ = A because of the orthorhombic distortion. 

The principal result conveyed in Figs. |2jb) and (c) is that 
the expected AFM order exists in the Ort-AFM phase. The 
question, however, is whether this AFM order is also associ- 
ated with the Ort-PM intermediate phase. Figures |2d) and 
(e) show [£ £ 0] scans through the (11 8)t charge and {\\ 
7)t magnetic peak positions at T = 133.3 K. Similar to what 
was found in our laboratory-based measurement [Fig. [Ha)] we 
observed four charge peaks [two outer peaks from Ort-AFM 
and two inner peaks from Ort-PM]. However, Fig. |2e) dis- 
plays only a single magnetic peak. The arrows in this panel 
denote the expected positions for magnetic peaks associated 
with each of the charge peaks in Fig.|2jd) and we see that the 
magnetic peak is found at a position that corresponds to one 
of the two outer peaks associated with the Ort-AFM phase. 
This allows us to conclude that the magnetic order is associ- 
ated only with the Ort-AFM phase. Taken together, the x-ray 
diffraction and XRMS measurements suggest that: (1) The 
orthorhombic distortion at Ts is best described as a second- 
order transition; (2) the structural and AFM transitions in the 
as-grown BaFe2As2 compound are separated in temperature 
by approximately 0.75 K and; (3) a first-order magnetic tran- 
sition at Tjy drives the discontinuity in the structural order 
parameter at 133.75 K. To further substantiate these conclu- 
sions, we now turn to a study of the evolution of the structural 
transition in electron-doped BaFe2 As2 compounds. 
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FIG. 2: (Color online) (a) The measured (11 8)t charge diffraction 
peak above the structural/magnetic transitions. Panels (b) and (c) 
show the (11 8)t charge peak and (\ \ 7)t magnetic peak at T 
= 130 K, well below the transition region. Panels (d) and (e) show 
the measured intensities at the (11 8)t charge peak and (55 7)t 
magnetic positions at T = 133.3 K. The arrows in (c) and (e) indicate 
the calculated magnetic peak positions corresponding to each of the 
charge peaks in (b) and (d), respectively. The fitted value for the 
width of the charge and magnetic peaks are the same. 
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FIG. 3: (Color online) (a) X-ray data, (b) resistance (black line) 
and its temperature derivative (blue line), (c) orthorhombic distor- 
tion and (d) FWHM of the split (11 10)t Bragg peaks measured for 
Ba(Fei_ 2 ,Co I )2As2 with x = 0.018. In panel (c) the structural and 
magnetic transition temperatures are marked. 



B. High-resolution x-ray diffraction and resistance 
measurements of Ba(Fei a; 00^)2 AS2 and Ba(Fei r Rh :E )2As2 

It has already been established that the substitution of Co 
or Rh for Fe in the parent BaFe2As2 compound results in 
new and interesting behavior.^ -12 ' 14-17 ' 23 ' 24 As doping is in- 
creased, both the structural and magnetic transitions are sup- 
pressed and split, with the structural transition occurring at 
higher temperature. In Ba(Fe 1 _ a .Co a; )2As2, for Co concen- 
trations 0.03 < x < 0.06, we enter a region of the phase di- 
agram where magnetism and superconductivity coexist and 



compete Within this region, the magnetic and structural 
transitions are well-separated in temperature, and continuous 
in nature (see for example, references I23I and |24|) . It is, there- 
fore, interesting to probe the behavior of the structural and 
AFM transitions, described above, as they evolve towards split 
second-order transitions with doping. 

To this end, we have performed high-resolution x-ray 
diffraction measurements on four doped BaFe2As2 sam- 
ples: Ba(Fei_ a; Co :E )2As2 for x = 0.018 and 0.047, and 
BafFei.^Rh^Asa for x = 0.012 and 0.040. Figures 
through [6] display the raw diffraction data, the orthorhombic 
distortion (<5) and diffraction peak widths derived from fits to 
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FIG. 4: (Color online) (a) X-ray data, (b) resistance (black line) 
and its temperature derivative (blue line), (c) orthorhombic distor- 
tion and (d) FWHM of the split (1 1 10)t Bragg peaks measured for 
Ba(Fei_ I Rh a: )2As2 with x = 0.012. In panel (c) the structural and 
magnetic transition temperatures are marked. 



FIG. 5: (Color online) (a) X-ray data, (b) resistance (black line) 
and its temperature derivative (blue line), (c) orthorhombic distor- 
tion and (d) FWHM of the split (1 1 10)t Bragg peaks measured for 
Ba(Fei_ 2 ,Co a: )2As2 with x = 0.047. In panel (c) the structural and 
magnetic transition temperatures are marked. 



the data, and the electrical resistance measured as a function 
of temperature. 

Turning first to the compounds at lower doping concentra- 
tions, Ba(Fe .982Coo.oi8)2As 2 and Ba(Fe .988Rho.oi2)2As 2 
(Figs. [3] and |4), below T$ both samples manifest a lattice 
distortion that evolves continuously as temperature is low- 
ered, until the onset of magnetic ordering where a step-like 
feature in the structural order parameter (5) is observed. At 
Tjy a distinct broadening of the split orthorhombic diffrac- 
tion peaks is evident over a narrow range in temperature. 
In contrast, the temperature dependence of the order param- 
eter and peak widths for the higher doping concentrations, 



Ba(Feo.953Co . 47) 2 As2 and Ba(Fe . g6 oRh .o4 ) 2 As2, is de- 
cidedly different near T/v (Figs. [5] and [6j. For these samples, 
the structural distortion evolves continuously, with only a mild 
kink in evidence at and without the attendant peak broad- 
ening at Tn- 

The differences between the lower and higher doping con- 
centrations are consistent with a change in the nature of the 
magnetic transition from first-order for low doping, to second- 
order for higher doping levels. At low-doping, as for the par- 
ent BaFe 2 As2 compound, there is a second-order transition 
from the tetragonal phase to the Ort-PM structure as temper- 
ature is decreased below Tg, The step in the orthorhombic 
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FIG. 6: (Color online) (a) X-ray data, (b) resistance (black line) 
and its temperature derivative (blue line), (c) orthorhombic distor- 
tion and (d) FWHM of the split (1 1 10)t Bragg peaks measured for 
Ba(Fei- I Rh a: )2As2 with x = 0.040. In panel (c) the structural and 
magnetic transition temperatures are marked. 



distortion (6) at Tjv is a consequence of the abrupt appearance 
of Ort-AFM coincident with a first-order AFM transition. We 
note that throughout this temperature range only two broad- 
ened peaks are observed in contrast to what was shown above 
for the parent compound. This is expected, however, since 
the larger separation of Ts and Tjv allows i5 to evolve to a 
value for the Ort-PM phase that is close to its magnitude for 
the Ort-AFM phase. The anomalous increase in the widths of 
the x-ray diffraction peaks at Tjv arises from the coexistence 
and near coincidence in position of the Ort-AFM and Ort-PM 
diffraction peaks over a narrow temperature range. For the 
higher doping levels, within our experimental resolution, the 
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FIG. 7: (Color online) Diagram showing the nature of the struc- 
tural and magnetic phase transitions for Ba(Fei- 2 ,Co a; )2As2 at Ts 
and Tjv, respectively. The thick line denotes a first-order transi- 
tion whereas the thinner lines represent second-order transitions. The 
crosses denote values for Ts and Tn determined from our measure- 
ments. The open circle denotes the approximate position of a tricrit- 
ical point as described in the text. 



absence of a distinct step in 6 or peak broadening at Tjv is 
consistent with second-order AFM and structural transitions 
as found for electron-doped BaFe2 As2 previously 

We summarize our results in Fig. [7] which displays a phase 
diagram for Ba(Fei_ ;E Co x )2As2 focusing on the concentra- 
tion range of the present study. The structural transition, over 
the entire range is second-order, whereas the magnetic tran- 
sition changes from first-order to second-order at a tricritical 
point as discussed below. 

IV. DISCUSSION 

To understand the existence, and estimate the location, 
of a magnetic tricritical point in the phase diagram of 
Ba (Fei_ x Co x ) 2 As2, we can, at first, rationalize the interplay 
between the magnetic and elastic degrees of freedom in terms 
of a simple Ginzburg-Landau model, similar to what was done 
by Cano et al. in Ref. |3J. We start from the effective free en- 
ergy: 

fas 2 U S r4\ , ( a m 2 . U m i\ x r 2 , n 

F eS =\—5 + —S ) + y — m +~ m J-Adm (1) 

with a s = a S fl (T - T s ,q), a m = a mfi (T - T Nfl ), and 
positive constants us, u m , A. Here, m is the antiferromag- 
netic order parameter and Ts,o, Tjv,o denote the bare struc- 
tural and magnetic transition temperatures respectively. For 
Ts.o < Tjffl, this model describes a simultaneous mag- 
netic/structural first-order transition. However, for Ts > 
Tjv.o> this model describes a second-order structural transi- 
tion at Ts = Ts o, followed by a magnetic transition at 



7 



0.0020 



0.0015 



0.0010 



< 0.0005 



0.0000 



"I r- 

\ 



1 



i 



Ba ( Fe i-x CO x) 2 AS 2 



_i I i_ 



0.00 0.01 0.02 0.03 0.04 0. 
Co concentration, x 



05 



FIG. 8: Jump of the orthorhombic order parameter AS = 
5ort-AFM — ^Ort-PM across the first-order magnetic transition, as 
function of x. The linear relationship AS oc (x — xtri) (dashed line) 
follows from the mean-field solution of Eq. [T] 
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FIG. 9: (Color online) The two magnetic ground states of the iron 
pnictides, characterized by stripes along two orthogonal directions, 
expressed in terms of two interpenetrating AFM sublattices with 
Neel vectors mi and rri2. Notice that the (x,y) coordinate system 
used here refers to the unit cell with two Fe atoms and is rotated 
by 45° with respect to the (a, b) coordinate system relative to the 
single-Fe atom unit cell. 



Tn (Tjsrfi < Tjy < Ts), which can be either first-order 
or second-order. Considering that Tn and Ts change with 
doping, the magnetic tricritical point takes place at the dop- 
ing concentration x% T i where u m agfi (Ts — Tn) = A 2 . Ex- 
perimentally, we know that (Tg — Tn) increases with dop- 
ing, x. Therefore, it is straightforward to conclude from the 
mean-field solution of Eq. (Q]) that, close to the magnetic tri- 
critical point, the jump in the orthorhombic order parameter 
across the first-order magnetic transition changes with doping 
as AS = c>ort-AFM - <5ort-PM oc (x - x tI i). Extrapolating 
this linear relation and using the values of AS from Figs. Q] 
and[3] we estimate that the magnetic tricritical point is located 
at Xt T i ~ 0.022, as shown in Fig. [8] 

The main issue with the model in Eq. (Q]l, however, is that it 
requires a fine tuning of the independent structural and mag- 
netic transition temperatures Ts,o and Tn.o across the phase 
diagram. In all of the phase diagrams of the iron pnictide com- 
pounds, it is observed that the two transition lines track each 
other very closely, even within the superconducting dome* 2 ^! 
This suggests that these two states are strongly coupled, rather 
than independent, as assumed by the previous model. To ad- 
dress this issue, it has been proposed that the particular mag- 
netic structure of the iron pnictides gives rise to emergent 
Ising-nematic degrees of freedom that couple to the lattice, 
inducing the tetragonal-to-orthorhombic phase transition^— 
In the magnetically ordered phase, there are two degenerate 
ground states characterized by in-plane spin stripes along each 
of the two orthogonal Fe-Fe bond directions. These ground 
states can be described in terms of two interpenetrating AFM 
sublattices with staggered magnetization mi and m.2, such 
that mi is either parallel or antiparallel to rri2 (see Fig. [9). 

Within this description, the magnetic free energy of the sys- 
tem F mag has contributions from each sublattice F{ and from 
the coupling between them, F%2. The former is given by: 



F '=2j J^~f X ^ (q) m% (q) ' m * ( ~ q) + 4 J ~ m * 

(2) 

where Xi (q) = Xo (^0 + ifo 2 — %T)z cos <7± c ^ is the static 

susceptibility of each sublattice, 34 u is a positive coupling 
constant, and the momentum q is measured relative to the 
magnetic ordering vector. Here, Xo 1 i s a magnetic energy 
scale, ro measures the distance to the magnetic critical point, 
a and c are the lattice parameters of the unit cell of volume v 
containing two Fe atoms, and r\ z is the inter-plane AFM cou- 
pling. The coupling between the two sublattices is given by: 




with g > 0, favoring configurations where mi and rri2 are 
either parallel or antiparallel. In a description of the magnet- 
ically ordered phase in terms of localized moments, this term 
originates from quantum and thermal spin fluctuations.— On 
the other hand, within an itinerant approach, where the mag- 
netic moments arise from the conduction electrons, the same 
term appears as a consequence of the ellipticity of the electron 
pockets^ 

The coupling (0) between the sublattices gives rise to an 
emergent Ising-nematic degree of freedom (p = mi ■ ni2, 37 
which may be finite even in the absence of magnetic order (i.e. 
(ip) 0, but (mj) = 0) as long as the magnetic fluctuations 
are strong enough^ 4 - A finite value, (mi • ma) ^ 0, breaks 
the Ising symmetry embedded in Eq. © and, consequently, 
the tetragonal symmetry. This can be seen explicitly through 
the magneto-elastic term: 



^mag-ci = A / S (mi • m 2 ) (4) 



where A > is the magneto-elastic coupling and 8 is the 
orthorhombic distortion. From the bilinear coupling of 5 
and if in Eq. both the nematic and structural transi- 

tions are simultaneous. This mechanism for the tetragonal-to- 
orthorhombic transition explains why the magnetic and struc- 
tural transitions track each other closely in all the phase di- 
agrams of the iron pnictides. Furthermore, it also explains 
several experimental observations, such as the softening of 
the lattice in the tetragonal phase and its hardening in the su- 
perconducting state, 34 as well as the suppression of the or- 
thorhombic distortion below the superconducting transition 
temperature M 

In the case where the elastic free energy is harmonic, F e i = 
C s S 2 /2 [where C s is the shear modulus] the only effect of 
the elastic degree of freedom is to renormalize the sublattice 
coupling constant g in Eq. ||3}, yielding g —> g + A 2 / C s ^ 
In a mean-field approach, this implies that the two transitions 
remain split and second order. Although fluctuations could 
induce a simultaneous first-order transition, 3 ^ it is unclear 
whether they could explain a second-order structural transi- 
tion split from a first-order magnetic transition, as our data for 
low doping levels in BaFe2As2 suggests. 

To account for our observations, we move beyond the har- 
monic lattice approximation to consider the effects of anhar- 
monic elastic terms (for more details on the formalism, see 
Ref. kll) . In the tetragonal phase, the most general form of the 
free energy can be written as F c \ = ^CijCiCj + ^Cijk€i£j€k, 
where Cy are given in terms of the six independent elastic 
constants and the strain components e< are: 

e\ = u xx + u yy + u zz ; e 4 = 2u yz 

£2= (u XX + Uyy - 2Uzz) ] ^ = 2u xz 

£3 = (uxx - Uyy) /\/2 ; e 6 = 2u xy (5) 

with Uij = (diUj + djUi) /2 and u = (u x ,u y ,u z ) denot- 
ing the displacement vector. In this notation, we have the 
orthorhombic distortion 6 = e^/2 and the shear modulus 
C s = 46*66 • Since we are interested in describing the transi- 
tion from the tetragonal to the orthorhombic phase, we retain 
only the essential anharmonic terms that contain e 6 : 

TjCWjj + 2 P 166£l ^266 e 2] ^6 + C45Q€4£ 5 €e (6) 

Minimization with respect to the other strain components 
yields an effective elastic free energy in terms only of ee = 26: 

F e} [6] = \c s 5 2 + -U5 A + lw6 6 (7) 
2 4 D 

where we included the sixth-order term W > to ensure sta- 
bility of the functional. Note that since: 41 

U_ _ T j _ (C^Cfgg - 2C , i 2 C l l66C , 266 + CuCl && ) 

16 2 (tf 11 q B -c? 2 ) 
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FIG. 10: (Color online) Phase diagram of the system with anhar- 
monic elastic terms. T denotes temperature, C s is the bare shear 
modulus, and Xa 1 is a magnetic energy scale. Thin (thick) lines re- 
fer to second-order (first-order) phase transitions, with the red (blue) 
lines denoting magnetic (structural) transitions; the simultaneous 
first-order transition line is denoted by the double line. We use the 
notation Ort' to emphasize that the orthorhombic distortion jumps 
across the first-order magnetic transition. The orange dotted line sig- 
nals the occurrence of a jump in both the magnetic and orthorhombic 
order parameters, without symmetry breaking. The open circle refers 
to the magnetic tricritical point, while the arrow indicates the value 
of C a for which we calculate the temperature dependence of both the 
magnetic and orthorhombic order parameter (see Fig. II 11 . 

the fourth-order coefficient can be negative, depending on the 
magnitudes of the anharmonic terms Cyjt. Here, Uq is the 
bare coefficient coming from higher-order quartic anharmonic 
terms. 

In what follows, we consider all elastic coefficients to be 
temperature independent, and that the only minimum of Eq. 
(0 is at 5 — 0. Thus, different from the model in Eq. ([TJ, the 
system has no intrinsic structural instability, and the elastic 
phase transition results solely from the magneto-elastic cou- 
pling in Eq. (0]). In the case of a harmonic lattice, it was shown 
that nematic fluctuations renormalize the shear modulus in the 
tetragonal phase, making it vanish when the magnetic corre- 
lation length achieves a threshold value. 34 Here, not only C s 
will be renormalized by nematic fluctuations, but also the an- 
harmonic term U in Eq. (0, which gives rise to a much richer 
phase diagram. 

To calculate the total free energy of the system, we use 
the 'mean-field 1/N approach' discussed elsewhere^i 3 ^ Ba- 
sically, we assume that the magnetic order parameter has 
N components and expand to leading order for large N. We 
then obtain self-consistent equations involving the magnetic 
correlation length £, the nematic order parameter f, the mag- 
netic order parameter m = |(mi)| = |(m2)|, and the or- 
thorhombic distortion 5. The latter is obtained by minimizing 
the effective elastic free energy F g = F c \ + F, where F is 
an implicit function of 8, arising from the 1/N solution of the 
magnetic problem with free energy F mag + F roa , g _ e i. Thus, F 
describes how magnetism changes the elastic free energy. 

We choose parameters that yield relative temperatures and 
jumps comparable to those observed experimentally (see Sec- 



tion IIIB)42 In particular, we take U < in Eq. ; once 
again, we stress that the values used for U do not introduce 
any local minimum other than 6 = in the bare F e i . To il- 
lustrate the richness of the resulting phase diagram, in Fig.fTOl 
we show the results obtained after fixing all parameters but 
the bare shear modulus C s - For smaller values of C s , the sys- 
tem undergoes a simultaneous first-order structural/magnetic 
transition from the tetragonal/paramagnetic phase to the or- 
thorhombic/antiferromagnetic phase. This corresponds to a 
direct first-order transition from the Tet-PM phase to the Ort- 
AFM phase, which has not been observed in the experiments 
reported here. 

As the bare shear modulus increases, the two transitions 
split: at higher temperatures, the system undergoes a second- 
order structural transition and then a. first-order magnetic tran- 
sition at lower temperature. The latter is accompanied by a 
discontinuity in the orthorhombic order parameter S due to 
the magneto-elastic coupling. This is precisely the sequence 
observed in our experiments described in Section III for the 
parent BaFe2As2 and doped compounds for low doping con- 
centrations (Tet-PM Ort-PM ->• Ort-AFM). Note that this is 
not another structural phase transition, but a consequence of 
the first-order character of the magnetic transition. To show 
this explicitly, in Fig. QT| we plot both S and m as function 
of temperature for the value of C s indicated by the arrow in 
Fig. [10] Not only is the relative size of the step comparable 
to that measured experimentally for BaFe2As2, but also the 
relative temperature at which the step occurs (see Fig. [T}b), 
where T s w 134.6 K and T N « 134 K). The discontinuity 
in the orthorhombic distortion accompanying the first-order 
magnetic transition is a very general feature that holds reg ard- 
less of the specific values of the parameters. Thus, it supports 
our interpretation that the experimental data in Fig. [JJon the 
parent compound, BaFe2As2, describe a second-order struc- 
tural transition followed by a first-order magnetic transition. 

Returning to the phase diagram of Fig. [10] we note that 
as the shear modulus is increased even further, the transi- 
tions remain split but the magnetic transition becomes second- 
order, as it is observed for higher doping concentrations in 
Ba(Fei_ 2; Co. c )2As2 and Ba(Fei_ a; Rh x )2As2. At low temper- 
atures, there is another line that marks a simultaneous jump 
in both the magnetic and the orthorhombic order parameter, 
without any symmetry breaking. However, we point out that 
in our 1/N approach we have not taken into account a key 
feature of the magnetically ordered state: the reconstruction 
of the Fermi surface. For instance, we note that x tl -i ~ 0.02 
is close to the composition where evidence for a Lifshitz tran- 
sition, below Tjy, in Ba(Fei_ a; Co K )2As2 has been reported 
by thermoelectric power, Hall coefficient measurements and 
angle-resolved photoemission. 43 ' 44 Therefore, features in our 
model deep in the magnetically ordered phase, such as this 
extended line, are likely to change once the reconstruction of 
the Fermi surface is considered. For instance, one possibil- 
ity is that this extended line terminates at a finite temperature 
critical point. 

Furthermore, although in the phase diagram of Fig. [TOlwe 
changed only the bare shear modulus, it is unlikely that this 
is the only modified parameter as doping is introduced in the 
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FIG. 11: (Color online) Magnetic (m, open symbols) and orthorhom- 
bic (5, filled symbols) order parameters as function of temperature T 
(in units of the structural transition temperature T s ) for the system 
indicated by the green arrow in the phase diagram of Fig. 1 101 



parent compound. In particular, the increase in the splitting 
between the transitions is much more modest in Fig. [TOlthan 
found experimentally (Fig. |7). It is possible, then, that other 
parameters controlling the splitting are also changed, such as 
the magneto-elastic coupling, A, and the inter-plane magnetic 
coupling, r\ z . The main objective of the phase diagram pre- 
sented here is to illustrate the various possible phase transi- 
tions once anharmonic elastic terms are taken into account. It 
is interesting to note that, in our simple phase diagram, sys- 
tems with softer lattices are more likely to display simultane- 
ous first-order transitions. Indeed, CaFe2As2, which is sig- 
nificantly softer than BaFe2As2, presents relatively strong si- 
multaneous first-order transitions.- 



V. SUMMARY 

Our high-resolution x-ray diffraction and XRMS studies of 
BaFe2As2 have provided several new insights concerning the 
nature of the structural and magnetic transitions in the fas- 
cinating 122 iron pnictide family. First, we find that the or- 
thorhombic distortion at Ts is best described as a second- 
order transition and that the structural and AFM transitions in 
the as-grown BaFe2As2 compound are separated in temper- 
ature by approximately 0.75 K. We propose that a first-order 
magnetic transition at Tjv drives the discontinuity in the struc- 
tural order parameter at 133.75 K, and this is consistent with 
our measurements of the evolution of the character of the tran- 
sitions in Co- and Rh-doped BaFe2As2. Using these results, 
we can provide an estimate of the position of a tricritical point 
in the phase diagram of Ba(Fei_ a; Co a ;)2As2. Finally, we em- 
ploy a mean-field approach to show that our measurements 
can be explained by the inclusion of an anharmonic term in the 
elastic free energy and magneto-elastic coupling in the form 
of an emergent Ising-nematic degree of freedom. 
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